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Abstract. We study a symmetric randomly moving line interacting by exclusion with a wall. 
We show that the expectation of the position of the line at the origin when it starts attached 
to the wall satisfies the following bounds: 

cit 1/4 <E&(0) < ^t 1 / 4 log t 

The result is obtained by comparison with a "free" process, a random line that has the same 
behavior but does not see the wall. The free process is isomorphic to the symmetric nearest 
neigbor one-dimensional simple exclusion process. The height at the origin in the interface 
model corresponds to the integrated flux of particles through the origin in the simple exclusion 
process. We compute explicitly the asymptotic variance of the flux and show that the probabil- 
ity that this flux exceeds Kt l / A logt is bounded above by const. t 2 ~ K . We have also performed 
numerical simulations, which indicate E£ t (0) 2 ~ t l ^ 2 \ogt as t — > oo. 
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1 Introduction 

We consider a process £ t on 

X = {£ G N z : - £(x + 1)| = 1, £(0) even} 

the space of trajectories of nearest neighbor random walks that stay non negative and such 
that at even "times" the walk visits even integers. 
The generator of the process is given by 

C f(0 = \ E ^ + A ^ 6 * * 0} [/(£ + A£(a:) 6 X ) - f(0] (1-1) 

X 

where 5 X is the infinite vector having 1 in the X th coordinate and zero on the others. The sum 
£ + ad x is understood coordinatewise. The discrete Laplacian A is defined by 

A£(x) := £(x + 1) - 2£(x) + £(z - 1). (1.2) 
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In words we can describe the dynamics as follows. The discrete Laplacian assumes only three 
values, —2, and 2. When the Laplacian is zero, the interface does not move. When it is —2 
or 2, at rate \ it makes a jump of length 2 in the same direction as the sign of the Laplacian. 
Over this motion we impose a restriction to keep the process in X: the interface cannot be 
negative, so we simply prohibit the jumps which violate the restriction. This is the meaning 
of the indicator function 1{£ + A£(x) 8 X > 0} in the generator. We can think the prohibition 
of becoming negative as the interaction by exclusion of the interface with a wall at —1. For 
shortness we call £t the wall process. 
Our main result is the following 

Theorem 1.3 Let £ t be the process with generator ( \1 . 1\ ) and initial flat configuration: 

f (x) := x(mod 2). (1.4) 

Then there exist positive constants C\ and c 2 such that 

cit 1/4 <E&(0) < c^ 4 log t (1.5) 

for sufficiently large t. 



Theorem |1 .3| catches the effect of the "entropic repulsion" in a stochastically moving interface 
interacting with a wall by exclusion. 

The line induced in R 2 by joining (x, £t(x)) to (x + 1, £t(x + 1)) for all x G Z has the same 
behavior as the interface between — l's and l's in a zero-temperature two-dimensional nearest- 
neighbors Ising model with a positive external field in the semiplane below the diagonal x = y 
and with initial condition "all plus" below the diagonal and "all minus" above it. See Section 
|^ for details. 

The equilibrium statistical mechanics of this model is well known. If one considers the 
generator C restricted to the box [—L,L] with boundary conditions £ t (— L) = £t(L) = 0, the 
invariant distribution is the uniform distribution in the set of nearest neighbors random walk 
trajectories starting at time — L at the origin, finishing at time L at the origin and being non 
negative for all intermediate times. Actually the uniform measure is even reversible for the 
process. But the uniform measure in this set corresponds to the law of a symmetric nearest 
neighbors random walk Xi conditioned to the set = Xl = 0, X,i > 0, i G [— L, L]}. Hence, 

the typical height of a configuration £ with the invariant law in the bulk of the box is 



a[rL]) ~ O(VrL) 

More precisely, the normalized process process (L 1 ' 2 £([Xr]), r G [—1,1]) converges as L — > oo 
to Brownian excursion on [—1, 1]; see Theorem 2.6 of Kaigh (1976). 

Many papers deal with the problem of entropic repulsion in Equilibrium Statistical Mechan- 
ics. The role of the entropic repulsion in the Gaussian free field was studied by Lebowitz and 
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Maes (1987), Bolthausen, Deuschel and Zeitouni (1995), Deuschel (1996) and Deuschel and 
Giacomin (1999). 

Entropic repulsion for Ising, SOS and related models was discussed in Bricmont, El Mel- 
louki and Frohlich (1986), Bricmont (1990), Holicky and Zahradmk (1993), Cesi and Mar- 
tinelli (1996), Lebowitz, Mazel and Suhov (1996), Dinaburg and Mazel (1994) and Ferrari and 
Martinez (1998). 

The exponent 1/4 for dynamic entropic repulsion was predicted by Lipowsky (1985) using 
scaling arguments. This exponent was then found numerically by Mon, Binder, Landau (1987), 
Albano, Binder, Heermann, Paul (1989-1992), see Binder (1990), De Coninck, Dunlop and 
Menu (1993). It has also been observed in real experiments by Bartelt, Goldberg, Einstein, 
Williams, Heyraud, Metois (1993). Further theoretical investigations of dynamics of lines, in 
relation to experiments can be found in Blagojevic, Duxbury (1999). 

Dynamic entropic repulsion for a line of finite extension L when t, L — > oo strongly depends 
on the ratio t/L 2 . The present paper deals with L = oo (analytical) or t/L 2 — > (numerical). 
The case t/L 2 = 0(1) has been studied by Funaki and Olla (2001). 

The exponent 1/4 also applies to the growth of fluctuations of an initially straight interface 
not interacting with the wall (see (p. 59|) below). For the Gaussian case, explicit computations 
were made by Abraham, Upton (1989), Abraham, Collet, De Coninck, Dunlop (1990). It was 
observed numerically in the two-dimensional Ising model by Stauffer, Landau (1989). 

The strategy to show Theorem [L3| is to compare the wall process £ t with a free process Q 
having the same local dynamics as £ t but not interacting with the wall. The free process lives 
in 

Xo = {(EZ z : \C(x) - ((x+ 1)1 = 1, C(0) = even} 

and its generator is 

4>/(C) = \ £[/(C + AC(s) 5 X ) - /(C)] (1.6) 

X 

In the next section we prove that with flat initial condition the variance of the height at the 
origin for the free process behaves as t x l 2 : 

lim r^VCtCO) = ±=. (1.7) 

We then couple the wall process and the free process in such a way that 

Ct(x) < £ t (x) (1.8) 

for all x and t. The free process has enough symmetry and, properly rescaled, has uniformly 
bounded in time exponential moments. With these, ( |1.7|) and (|1.8| ), we get the lower bound in 

& 

The idea for the upperbound is to consider a family of free processes with initial condition 
depending on t: 

(% t (x) = ( (x) + a t 
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(a flat interface of height at). Then we fix a t = ct 1 ^ 4 logt, the constant c to be determined later 
and exhibit a coupling under which 

e.(o) <c(°) (i-9) 

for all s < t with large probability. Combined with (|1.7| ), inequality ( |1.9|) is the key for the 
upperbound in ( |1.5| ). The existence of exponential moments (mentioned above) yields the 
moderate deviations result needed here. 

The control of the fluctuations of the position at the origin of the free process is obtained by 
an isomorphism between the free process and the one-dimensional symmetric nearest-neighbor 
simple exclusion process. Under this map, Ct(0) — 2J t , where J t is the integrated flux of 
particles at the origin for the exclusion process. We compute explicitly the asymptotic variance 
of the integrated flux for the flat initial condition in Theorem |2.17| below and obtain 

, VJ t 1 , 
lim —± = —=. 1.10 

De Masi and Ferrari (1985) proved that the asymptotic variance of the integrated flux when 
the initial configuration is distributed according to a product measure with density 1/2 is given 
by 

li m ^ = i (1 .n) 

y/t 2V2TT 



which is strictly bigger than (|1 . 10 ) . When the initial density is p, the asymptotic variance 



is given by p(l — p)^2/ir. The method to show ( 1.101) and ( 1.11|) is based on duality and 



comparison with systems of independent particles and it is inspired by Arratia (1983), who 
used these tools to compute the variance of a tagged particle for the process starting with an 
(invariant) product measure. However a modification of Arratia's proof is needed in ( l.lOj) due 
to the deterministic character of the initial configuration. 

The study of the flux in the simple exclusion process is done in Section 2. In Section 3 
we prove Theorem |1.3| . Section 4 is devoted to numerical simulations. Section 5 shows the 
equivalence to the dynamics of a particular zero temperature Ising model interface. 

2 The free process and simple exclusion 

The simple exclusion process lives in {0, 1} Z and its generator is 



where 



r)(y) if y ^ x, x + 1 
r] x ' x+1 {y) = { r}{x + 1) ify = x (2.2) 
Tj{x) if y — x + 1 
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It is convenient to construct the processes using the Harris graphical construction. 

Harris graphical construction Let (N t (x) : x G Z) be a family of independent Poisson 
processes of rate |. For each x, N t (x) counts the number of Poisson events associated to x in 
the time interval [0,t]. Denote dN t (x) = ljthere is a Poisson event associated to x at time 
t] = lim/j^ (A^ t (x) — N t _ h (x)). Let r\ t be the process defined by 

d7] t (x) = (r] t (x - 1) - 7] t (x))dN t (x) + (r] t (x + 1) - T] t (x))dN t (x + 1). (2.3) 

The process is well defined because for each finite time t the value of the process in a finite box 
can be determined by looking at only a finite but random number of Poisson events and initial 
values. See for instance Arratia (1983). In words, the motion can be described as follows. The 
Poisson marks of N t (x) are associated to the bond (x,x + 1) and each time a Poisson mark 
occurs, the contents of the associated bond are interchanged. It is immediate to show that this 
process has generator ( |2.1| ). 

Stirring particles. To introduce the notion of duality and to deal with the flux of particles 
it is convenient to follow the "stirring particles" as defined by Arratia (1983). Let X* be the 
position at time t determined by Xq = x and the equations 

dX x t = dN t (Xf) - dN t {X? - 1) (2.4) 

So that, each time a Poisson mark associated to one of the neighboring bonds of a particle 
occurs, the particle jumps across the bond. Of course, if both extremes of a bond are occupied, 
the particles jump simultaneously, respecting the exclusion condition "at most one particle per 
site" . For each t > the (random) map 

x ^ X x t (2.5) 

is a bijection of Z in Z. The (marginal) law of X* is a symmetric nearest neighbor random 
walk starting at x. 

Duality. Let y i— > D\ be the inverse map defined by x = D\ if and only if y = Xf. The 
following "duality formula" holds immediately 

Vt (y) = 77o(A a ) (2-6) 

So, 

Hih{v) = Ilva(Df) (2.7) 

y&A yeA 

Notice that for a finite set of sites A, {D% : y e A} has the same one-time marginal as a simple 
exclusion process with initial condition A (here we are identifying the configuration rj with the 
set {x : f](x) = 1}). When A = {y} (contains only one site), the one-time marginal D\ has the 
same law as X\ for all t > 0. 
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Integration by parts formula Consider (Y t \ Zf) independent random walks with the same 
marginals as the stirring process (JQ, Xf). The generator of the process (Y t \ Z ] t ) is the following: 



Uf(i,j) = l E [f((i,j) + e)-f(i,j)} (2.1 



2 

eGZ 2 :|e|=l 

and the generator of the process X t is 

(\Y.e^:\e\=Af^3)+e) ~ f(i,j)] if i - j > 1 

^J)=< i/(j,i) + |/(i + l,i) + |/(z,i-l)-|/(i,i) ifi-j = -l (2.9) 
I I/O", + - 1, J) + lf(i,3 + 1) - J) if i -j = 1 

Hence, for i ^ j, 

J) - V7M = " ~j\ = !}(/(*, J) + " /M - /(J, J)) (2-10) 
Let £/f and Vt be the semigroups generated by U and respectively. Let / : Z 2 — > R. Then 

E/(x t i ,x/)-E/(y^,z|) = [y t -i7 t ]/(i > j) = fv a [u-v]u t _ s f{i,j)ds (2.11) 

./o 

where the last identity is the integration by parts formula (see Liggett (1985) Proposition 8.1.7). 
Now, using (|2.10|) to compute ( |2.11|) we get for i ^ j: 

Ef(Xi,Xi)-Ef(Yt,Zi) 

' 1 dsE(l{\Xi-Xt\ = l} (2.12) 



x [f(Y t xt -%z?*-°) + m xi -%z? L °) - f(Y t xi -%zp-°) - m xL %z?*-°)] 

This identity will be used in the sequel. 

Flux. Let J t be the integrated flux of rj particles through the point —1/2 in the exclusion 
process: 

J t := > 0} - J2vo(x)l{X? < 0} (2.13) 

x<0 x>0 

where Xf is the position at time t of the exclusion particle that at time zero was at position x. 
Replacing ( |1.4p in ( |2.13| ), we write 

Jt := E W > °) " E < °>- ( 2 - 14 ) 

Jt is almost symmetric. Let 

^ := l{Xf > 0}; H' t := £ !«* < "!}; * == E W < °>- ( 2 ' 15 ) 

j<0 i>0 i>0 
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Then, clearly, 

J t = H t -I t ; H t ~H' t - \H' t -I t \<l, (2.16) 

where ~ means identity in distribution and is justified in this case by spatial and distributional 
symmetry. 



Theorem 2.17 For the simple exclusion process with generator (\2.1\ ) and initial condition 
r/ = 1 — 2£ , as defined in t \l-4l ), 

VJ < 1 (2.18) 



lim _ — _ . 



Proof. Working from (|2.14j) , we get 

E(J t ) 2 = J^P(Xf > 0) + ^P(Xf < 0) + 2 E F ( X ? - °' X t j - °) 

i<0 i>0 i<0 i<jr'<0 

+ 2 E p ( x21 < °' x ? < °) ~ 2 E E p ( x2i - °' x ? < °) ( 2 - 19 ) 

i>0 i>j>0 i<0 j>0 

(EJ t ) 2 = ^P 2 (Xf > 0) + J> 2 (X 2i < 0) + 2 E ¥ ( X ? - °) ¥ ( X t 3 > 0) 

i<0 t>0 i<0 i<j<0 

+ p ( x2i < °) p ( x2i < °) ~ 2 E E p ( x * 2i - °) p ( x * 2i < °) ( 2 - 2 °) 



*>0 i>j>0 

Immediately we have: 



j<0 j>0 



nxf > o) + J2 nxr <o) = j2 n*i > o) 



(2.21) 



j<0 i>0 i<0 

and analogously for the P 2 terms in ( p.20| ). Using 

F(AB)-F(A)¥(B) = -(P(A B c ) - F(A) ¥(B C )) = F(A C B c ) - F(A C ) ¥(B C 
we get 



where 

and 
St 



YJ t = V t + S t , 

v t = j2nxi>o)-J2 F2 ( x t>°) 



(2.22) 

(2.23) 
(2.24) 



i<0 



i<0 



E + E + 2 E ) ( p ( x2t ^ °« x * 2j ^ °) - p ( x21 ^ °) P ^ J ' ^ °) 

i,j<0,i^j i,j>0,i^j *<0,j>0 



J^(P(Xf > 0, X 2j > 0) - P(Xf > 0) P(X 2j > 0)) 

i+3 



(2.25) 



Since P(X| > 0)P(X/ > 0) = P(l? > 0, Z\ > 0), we can use (gig) with = l{i > 

0,j>0} toget 

P(X| > 0, X/ > 0) - F(Xl > 0) P(X/ > 0) 
= ~\ jf E P (^' = ^'^+ !}) ( P (*£. ^ °) -W-t 1 > 0)) 2 rf S (2.26) 

See also Theorem 2 of Ferrari, Galves and Landim (2000) for a probabilistic proof of the previous 
identity. Translation invariance and self-duality of (X*, X%) implies that ( |2.26| ) equals 



- \ f E p « x - = {< - j - y» p W-. = 

From ( p725D and (ggg)-(gg^), 



(2.27) 



ft 



Since X°±Xl 



\ f E p2 (^°- = f ) E p « x - x ^ = { 2i - v. - ?/}) 



(2.28) 



E P (W> X s 1 } = {2z-y,2j- 2/ }) = 2P(X s 1 (mod2) = X°(mod2) = y(mod2)) (2.29) 

i¥=3 

Let A y = {(z, j) G Z 2 : i(mod2) = j(mod2) = y(mod2)}. We show below that 

lim P((X S °,X!) e A y ) = 1/4. (2.30) 



uniformly in y. Hence, 



lim t-^St = -\ lim r 1 / 2 f Y P 2 (F S ° = y) ds 



(2.3i; 



Let Z° be an independent copy of Y t °. Since Y^y^^s = v) = F ( Y t° ~ Z t = °)> changing 
variables the above limit equals 

lim l\stfl^t ~ Z°st = 0) ^ = f lim ( S t) 1/2 P(F s ° t - Z° st = 0) ^, (2.32) 

where the interchange of the limit and the integral are guaranteed by the local central limit 
theorem for (Y® — Zf), which is a symmetric random walk of rate 2. This also implies that 
( 2.32|) equals 



ds 
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(2.33) 



We conclude that 



lim t~ 1/2 S t = =. (2.34) 



To compute V t notice that 

J2nXl>0) = i) = E((X°)+); (2.35) 

^p 2 (x;>o) = ]Tp 2 (x t °>*) = ^P(y t °AZ ( °>i) = E[(y t °AZ f °) + ], (2.36) 

i<0 i>0 i>0 

Thus 

limr^H = E(I + )-EpAl') + ] = — ^, (2.37) 

where X and X' are i.i.d. standard normals. 

Finally , substituting ( g3§ and (g57| ) in ( ggg) we get (gTjg ). □ 

Proof of ( |2.30 ) The continuous time Markov chain (3^°(mod2), F t 1 (mod2)) converges expo- 
nentially fast to the uniform distribution in {(0, 1), (1, 0), (0, 0), (1, 1)}. This implies that there 
exist positive constants C\, Ci such that 

|P((F S \ Z{) G A y ) - 1/4| < C ie - C2t . (2.38) 

uniformly in i,j,y. Writing f y (i,j) := l{(i,j) G A y } and using (|2.12|) we get 

\V((X°,X})e A y )-F((Y t °,z})eA y )\ 
< ±J*dsE(l{\X s -Xl\ = l} (2.39) 

x [f y (Y t X !-%Z^-°) + f y {Y^%Z? L °) - f y (Y t X? -%zf-°) - f y {Y t Xl -%zf^)] 



< 



2 f ds ¥{\X\ - X°| = 1) d e - C2(t ' s) (2.40) 
Jo 

(using ( |2.38|) ). Now \X} — X°| is a Markov chain in {1,2, . . .} with rates p(l,2) = p(x,x + 
1) = p(x,x — 1) = 1/2, x > 1. It can be easily coupled to a a Markov chain in {0, 1,2, . . .} 
starting in 0, say R®, with rates p(0, 1) = 1, p(x,x + 1) = p(x, x — 1) = 1/2, x > in such 
a way that |Xj — X°| > R® for all s. Since R® is a simple symmetric random walk reflected 
at the origin, we get that lim^^ P(|X| - X°| = 1) < lim^^ P(|i?°| < 1) = and thus, 
from ( |2.40| ) and dominated convergence (after a change of variables s — > t — s'), it follows that 
lim^ oo P((X s ,X s 1 ) G A y ) = hm^ oo P((y s ,Z s 1 ) G A y ) = 1/4. □ 

Lemma 2.41 Let H t be as in ( ^Wj ) and H t = t^/ 4 [H t - E(H t )\. Then for allXeR 

limsupE(e A ^) < e sA2/2 , (2.42) 

t^oo 

where s = l/^/2n. 
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Proof. 



E(e 



E[exp {\t-V*Y, i<0 l{X* >0})] 



^2.43) 



exp (At-V^. <o P(xf >0)) 

We will show that the quotient in ( 2.43|) is bounded above by a constant. For that, we need 
to evaluate the expected value in that equation. Let A > 0. We will argue below that 



E exp f Ar 1/4 ^l{Xf > 0} ) 

The last expectation equal; 
1 + [exp(At- 1/4 ) - 



< nE[exp(^ 1/4l {^>0})]- 

i<0 



The last expectation equals 

1 + [exp(At- 1/4 ) - l]P(Xf > 0) = 1 + [Ar 1/4 + (A 2 r 1/2 /2) + o^ 1 / 2 )]?^ > 0), 
for all t large enough. The last expression is bounded above by 

exp {[At- 1/4 + (AV 1/2 /2) + o(r 1/2 )}F(X 2i > 0)} . 
Substituting into the right hand side of ( |2.44j ), we get 

exp <! [(AV 1/2 /2) + o(r 1/2 )] J^P(Xf > 0) 



(2.44) 

(2.45) 
(2.46) 

(2.47) 



i<0 



as an upper bound for the quotient in (|2.43|) . It is not difficult to see that the expression on 
the exponent in (|2.47|) converges to e sX I 2 . 

To finish the argument for A > 0, we have to justify the inequality ( |2.44| ). That follows from 
taking limits as M — ► — oo (and using monotone convergence) on the respective inequalities 
gotten by replacing the infinite sums by ^A/<i<o- These are justified by the fact that the 
functions exp(t -1 / 4 ^ M<i<0 1{X 2 * > 0}) are bounded, symmetric and positive definite for all 
M < 0. The inequalities then follow from Proposition 1.7, Chapter VIII of Liggett (1985). 

For the case A < 0, we use the identity 



>0} >0) 

i<0 i<0 

and a similar argument as above. □ 



£l{X*<0}-$>(X*<0) 



i<0 



i<0 



f2.48) 



Lemma 2.49 Let t]q be given by the flat condition as in Theorem \2. 1 % Then 

.E( e l J */* 1/4 l)<oo. 



sup . 

t>0 



Furthermore for all K > and all t large enough 

F(\J t \> Kt 1/4 hgt) < cr K , 

where c is a constant. 



(2.50) 



(2.51) 
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Proof. The bound (|2.50|) follows straightforwardly from Lemma (|2.41|) and relations ( 2.16 ). 

From the relations ( |2.16| ), to show ( |2.51| ) it is enough to prove the result with \H t — ~E(H t )\ 
replacing \J t \ in (|2.51| ) (the constant of course does not need to be the same). We have 

F(\H t \ > Kt^Hogt) = F(\H t \ > \ogt K ) < c't~ K (2.52) 

where the last inequality follows from the exponential Markov inequality and d = sup <>0 E(e' i/i ') 
is finite by (ggp . □ 

Graphical construction of free process. Let be the process defined by 

d( t (x) = A( t (x)dN t (x), (2.53) 



where the discrete Laplacian A was defined in ( |1.2|) . In words, each time a Poisson mark of the 
process N t (x) occurs, the height at x at time t decreases or increases two unities, according to 
the value of the Laplacian at this point at this time; if the Laplacian vanishes, no jump occurs. 
This process has generator (|1.6|). 

Lemma 2.54 Let r] (x) = (o(x + 1) — (o(x). Then 

Vt (x) = ( t (x + l)-( t (x) (2.55) 

where the processes ( t o,nd r\ t are defined by ft273j) and ^2. 5$ ) and have initial conditions r] and 
(o respectively. Furthermore, 

0(0) - Co(0) = 2J t . (2.56) 

Proof. Notice that from ( |2.55| ), 

AC t (0) = 2( Vt (-l) - Tk(0)) (2.57) 

Assume that there is a mark of the process N t (— 1) at time t. Then ( |2.57|) , ( p.3|) and (|2.53|) 
imply that if rj t (— 1) — rj t (0) = no changes occur neither for rj t (— 1), rj t (0) nor for Ct(0); if 
i] t (—l) — rj t (0) = 1, an exclusion particle jumps from —1 to and the free process at the origin 
jumps two units up; if rjt(— 1) — T)t(0) = —1, an exclusion particle jumps from to —1 and the 
free process at the origin jumps two units down. Identity fl2.56|) follows from ( |2.55|) . □ 



Lemma 2.58 Let Q be the free process with flat initial condition Then 

lim ^# = -L (2.59) 

supE( e lfi(0)/tl/41 ) < oo (2.60) 

and for all K > and all t large enough 

P(|Ct(0)| > i^t 1/4 logt) < ct~ K . (2.61) 
Proof. It follows from identity ( |2.56| ), the limit ( |2. 18| ) and the bounds ( |2.50| ) and ( p.51| ). □ 
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3 Coupling the wall and the free processes 



We construct graphically the wall process which simultaneously provides another graphical 
construction for the free process. Under this construction the wall process dominates the free 
one. We consider two independent families of Poisson processes with the same law as N t (x) 
called Nf(x) and Nf(x), to be used for upwards and downwards jumps, respectively. The 
process satisfying the equations 

d£ t {x) = A£ t {x)l{A£ t {x)>0}dN+(x) + A£ t (x) l{A£ t (x) < 0, £ t (x) + A£ t (x) > 0}dNf(x) 

(3.1) 

has generator (|1.1|) . The process ( t satisfying 

d( t (x) = A&(z) l{A£ t (x) > 0} dN+(x) + A&(z) l{A&(x) < 0} dN t ~(x). (3.2) 
has generator (|1.6|). 

In words, when a time event of the process N^~(x) occurs at time t, the process £t at site x 
and time t jumps two units upwards if A^(x) > 0. When a time event of the process Nf(x) 
occurs at time t, the process £t at site x and time t jumps two units downwards if A^ t (x) < and 
the wall condition £ t (x) + A£ t (x) > holds. The process satisfying (|3.2| ) follows the same marks 
in the same manner but ignoring the wall condition. The difference with the process satisfying 
(|2.53|) is that in this case the Poisson events N t are used for both upwards and downwards 



jumps; this construction is not attractive in the sense that it does not satisfy ( |3.4| ) below. 

Let r be a non negative integer and and Q be the processes defined by (|3.1|) and ( |3.2| ) 
but with initial condition 

£r( x ) = Q( x ) = r + x(mod2) (3.3) 

Notice that Ct° and £t as defined in ( |2.53| ) have the same law but are different processes. The 
processes £ t and defined by ( [3.1|) and the same initial condition satisfy 

6(0) < Ct(0) (3.4) 

for all r > 0. This joint construction corresponds to what Liggett (1985) calls basic coupling. 

Lemma 3.5 There exists a constant c > such that for any K > and t > 

P(&(0) > 2Kt 1/4 log t) < ct 2 ~ K (3.6) 

Proof. Let a t = 2Kt 1 ^ 4 \ogt. Take r > and write 

P(&(0)>a t ) < P(^(0)>a t ) 

= P(e[(0) > a t , Ct(0) = C(0)) + P(£[(0) > a t , £(0) + (0)) 

< P(C(0) > a t ) + F(C t (0) ? C(0)) (3.7) 

The first term in (|3.7|) will be bounded using Corollary |2.58| . To bound the second term notice 
that if the interacting process and the free process differ at the origin this is due to a collision 
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of the interacting process with the wall at some point x that separate the two processes at x at 
some time s; the discrepancy then propagates and arrives to zero by time t. We fix an a > 
and separate the discrepancies in two classes: those that come from the interval [—at, at] and 
those that come from outside this interval. If in the time interval [0, t] the free process does not 
touch the wall in the space interval [—at, at], then the discrepancy must come from outside. 
Hence, 

{£[(0) ^ Ct (0)} C {(Z(x) < for some s E [0,t], x e [-at, at]} 

U { a discrepancy from [—at, at] c reaches up to time t} (3.8) 

Observe that the law of Q(x) — r is the same as the law of Ct(0) an d that P(£ s (0) < — r) < 
(0) > r) due to the initial condition being non-negative. Hence, fixing 



r = a t /2, (3.9) 

the probability of the first event in the right hand side of ( [3.8|) is bounded by 

(2at + 1) P(C(0) > at/2 for some s e [0,t]) (3.10) 

From (|2.6CI|) and the exponential Markov inequality, we have that sup s<4 P(£ s (0) > r) < ct~ K 
for some constant c, so we can bound ( |3.10| ) with 

(2a* + 1) I ct- K ds < (2a + l)ct 2 - K (3.11) 
Jo 



To bound the probability of the second event in the right hand side of ( |3.8| ) notice that 
discrepancies cannot travel faster than N t , a Poisson process of parameter 1. Hence 

P(a discrepancy from [-at,at] c reaches up to time *) < 2F(N t > at) < 2e~ t(a+1 ~ e) (3.12) 

using the exponential Chebyshev inequality. Fixing a = 2, and using the bounds ( p.lOj ) and 
( |3 • 1 2| ) , the probability of ( |3.8j ) is bounded by 

4ct 2 - K + 2e"' (3 - e) < c't 2 - K (3.13) 

for some constant d and sufficiently large t. □ 

Proof of Theorem |1.3| . It follows straightforwardly from ( |2.60| ) that C, 2 is uniformly inte- 
grable, where (t = ^ _1//4 Ct(0)- This, together with (|2.59|) , implies the lower bound in ( |1.5|) , as 
we will see now. Indeed, (Q implies that r^E^O) > E|C*|. Now, 

YQ < EC, 2 = E(( 2 ; C 2 < M 2 ) + E(C 2 ; ( 2 > M 2 ) < ME\( t \ + e M (3.14) 

uniformly in t, where M is an arbitrary positive number and €m —* as M — > oo. Thus 
t-^E^tiO) > (V( t - e M )/M. We conclude that 



liminf t~ 1/4 E&(0) > sup ( -= - e M ) /M > 0. 

M>0 \ V n J 



(3.15) 
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For the upperbound, we use [2.61| to obtain 



= E P ^(°) > kt^hgt) <4 + J2c't 2 - k/2 <c 2 < oo (3.16) 

°^ fc>0 fc>5 

for some constant C2 < 00. □ 

4 Numerical simulation 

We have simulated the processes ( t and £ t numerically, using various pseudo-random number 
generators. The interface is of length L with periodic boundary conditions, so that the processes 
live in 

Xol = {( e Z Z / LZ : |C(x) - ((x + 1)1 = 1, C(0) even} (4.1) 

or 

X L = {£ E N Z/LZ : \£(x) -£(x + l)\ = 1, £(0) even} (4.2) 

Time is an integer multiple of 2L~ l , i.e. t G (2L^ 1 )N. For each time step, a site x is chosen ran- 
domly according to the uniform measure on Z/LZ, and the interface is updated with the same 
rules as in the continuous time version of the processes. The transition operator corresponding 
to a time step 5t = 2L~ l is 

T f(C) = /(C) + L- 1 [f(C + AC(*) S x ) - /(C)] (4.3) 



or 



Tf(0 = f(0 + L- 1 ]T l{e + A^(x) 5 X > 0} [/(C + Af(x) 5,) - /(0] (4-4) 



Tf(£) is the expected value of the function / evaluated at time 2L _1 (after one step) when the 
initial configuration is £ for the discretized version of the process. The same interpretation is 
valid for T /(C). By abuse of notation we call the discrete time versions of the process £ f and 
( t as we did for the continuous time versions. Notice that 

^m^im =cm h whM =w(0 , (4 . 5) 

The processes £t and Ct are coupled in the simplest possible way: the same random sequence of 
sites are used for both. Notice however that this coupling is different from the one described in 
Section 3 (in particular it is not attractive in the sense that it does not necessarily satisfy ( |1.8|) 
but it is faster). For L finite the discrete time and continuous time processes can be identified 
up to a time change, using the ordered sequence of updated sites. The random time involved 
in the time change has fluctuations which should be negligible for our purposes. 
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The numerical samples for the data shown below were drawn using either the Mersenne 
Twister pseudorandom integer generator, see Matsumoto, Nishimura (1998), or the R250 pseu- 
dorandom generator, see Kirkpatrick, Stoll (1981). The length L is 10 6 or 2 20 = 1024 2 and 
time runs up to 2.10 6 or 2 21 . The number of calls to the generator for the realization of one 
sample of length L up to time t is L.t/2 < 10 12 , which of course is much less than the period of 
the generator (a necessary but not sufficient condition for reliability). We compute empirical 
averages 

e t =L- iy £^)\ (? = L-iJ2(t(x) 2 (4-6) 

X X 

and empirical distribution functions 

/,(n) = L- 1 ^l{6(^)=n}, f 0it (n) = L- l Y / nCt(x)=n} , neZ (4.7) 

X X 

scaled into 

Ms) = t'/'L- 1 £ l{t-V%{x) = s} , Ms) = t^L- 1 l{t- l/ %(x) = s}, se t-^Z 

X X 

(4.8) 

which, extended to s G R approximate the Schwartz distributions 

Us) = L- 1 S(t- lf %(x) - s) , Ms) = L- 1 ty-V'Ctix) -s), seR (4.9) 

X X 

where L _1 ^ x is an ersatz for the expectation over a real random variable, limit of i -1 / 4 £ t (x). 

The processes were studied for time t < 2L, whereas the effect of finite size with periodic 
boundary conditions is expected to be visible only after a time of order L 2 , the relaxation 
time of an interface of length L. The law of large numbers in empirical averages as above is 
believed to be at work with an effective number of weakly dependent variables of order L/t 1 / 2 : 
the interface at time t can be thought of as a collection of Ljt 1 ! 2 segments of length t 1//2 , the 
different segments being weakly dependent. For t — L and one sample of the process, we have 
only t l l 2 independent segments, hence an expected relative statistical error of order t" 1 / 4 . This 
explains the more erratic behaviour at larger times in Fig. 1. 



15 



2 



1.5 



0.5 





'' \ 1 t (averaged) 






MT 

JXiJU 

fit 


£ / * 1/4 

C^/ I (maximum probability) 




^ ^-V/-V-X- r -''^ --• ■- ...v -.:V^.>W^/ 

-1: 

f 

1 
1 
1 




1 
1 




C / t (averaged) 
i i i i 



500000 1e+06 1.5e+06 2e+06 

Fig. 1: from top to bottom, as function of time : t^ 1 / 2 ^ 2 together with best fit 1.62 + 0.024 logi ; the value 
of s where 4>t(s) is maximum, together with best fit -y/0.55 + 0.057 logt ; and t~ x l 2 Q together with the exact 
asymptotic value a/1/7t. Graphs labelled MT are averages over 6 runs with the MT random generator with 
different seeds. Graphs labelled R250 are averages over 5 runs with the R250 random generator with different 
seeds. Interface length is L — 2 20 or L = 10 6 . 

The numerical experiment clearly favors an asymtotic behavior E£ t (0) 2 ~ t 1 / 2 log t as t — > oo. 
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Fig. 2: <j)t(s) and 0o,t(s) at large times, from the same data as in Fig. 1. 

Figure 2 shows the scaled empirical distribution functions at various large times. Clearly 
t~ l ^Ct{ x ) converges to a centered Gaussian random variable as expected. The distributon 
function of t~ l l A ^ t {x) is markedly asymmetrical. Zooming around s = indicates 0oo(O) = 
0^(0) =0 and C(0)>0, and 

/,(0) = L- 1 ^l{6(x)=0}~r 1 /2 (4.10) 

X 

5 Interface of the Ising model at zero temperature 

In this section we explain the relation of our model with the interface of a particular Ising 
model at zero temperature. Let the "inverse temperature" (3 > and a t G { — l,+l} z2 be the 
Ising model with generator 

= ^£ [/(**)-/(*)] 
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with o x {z) — a(z) for z ^ x G Z 2 , cr x (x) = — a(x) and Cp(x, a) are the Glauber rates 



e -PH(o*) 

with Hamiltonian 

for some positive magnetic field h > 0. Consider the case /3 = oo and assume that the starting 
configuration <jq is "all ones" below the diagonal and "all minus ones" above or in the diagonal: 

={!l; (5 - 1) 

In this case for all t the configuration 07 has the property that all sites have either exactly two 
or no neighbor with opposite sign; furthermore, only sites above or in the diagonal may be 
negative. As a consequence, for all t > the rates c(x, a t ) are positive only for sites x above 
or in the diagonal for which there are exacly two neighboring sites with different sign: under 
initial condition (15.11), 



( x? a< ) = «) 1/2 if E y; | r x|=i H^tiv) + 0tO)} = 2 and x i < x 2 
otherwise 



To get the wall process of Theorem [T7| from the above dynamics with initial condition (p.l|) , 
we first rotate the lattice by —45° and multiply by a/2, that is, we perform the transformation 
R : Z 2 — » R(x,y) = (x + y,x — y), where Z% := {(x, y) G Z 2 : x + y is even} is the 
even sublattice of Z 2 . The above dynamics then induces a dynamics in { — 1,+1} Z 2 given by 
&t(z) = cr t (i? _1 z), z G Z 2 ,t > 0. Defining £t(ar) := min{y : G Z 2 and a t (x,y) = — 1}, 

x G Z, t > 0, we have that the wall process £.(•) with generator ( |1 . 1| ) and initial configuration 
( |1.4| ) has the same law as £.(•). 
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